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By Donald Bogart and Michael Е. табо Я 
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SUMMARY 


Results are presented of two-group criticality caYculetiené made for 
sodium-hydroxide-cooled, moderated, and reflected геабфогв ғов various 
concentrations and compositions of reactor fuel- element strüctural 
material. These specific criticality results are presenti in & general- 
ized manner to permit rapid evaluation of the criticality’ requirements for 
a wide range of other structural material compositiongsand concentrations 
which may be of interest from considerations of ыа нау Strength, апа 
heat-transfer-surfa&ce requirements. 


Based on turbojet-engihe cycle operating conditions optimized for 
minimum airplane gross weight, the maximum reactor fuel-element and 
coolant temperatures are related to the reactor heat release and airplane 
gross weight for а range of the reactor heat-transfer variables for flight 
at altitudes of 30,000 and 50,000 feet and Mach numbers of 0.9 and 1.5. 
For the calculations, airplane lift-drag ratio is assumed constant at 6.5 
for supersonic and 18 for subsonic flight. The weight of shield pius 
reactor plus payload plus auxiliaries (herein designated as Wk) is assumed 
constant at two different values namely, 100,000 and 150,000 pounds; the 
calculated results for 100,000 pounds are considered representative for the 
divided-type shadow shield, the results for 150,000 pounds apply for the 
unit bulk shield. 


The results provide a basis for compromise of the advantages of higher 
cyele efficiencies (and hence lower airplane gross weights, reactor heat 
releases, and engine air flows) attainable at high reactor fuel-element 
temperatures with the advantages of higher reactor-materlal strength and 
corrosion resistance attainable at Low fuel-element temperatures. 


For subsonic flight at altitudes of 30,000 and 50,000 feet and for 
values of Wg of 100,000 and 150,000 pounds, maximum reactor fuel-element 
temperatures of the order 11009 to 12009 F сап be maintained with а reactor 
core diameter of 2 feet. The airplane gross weight and reactor heat 
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release are approximately 200,000 pounds and 70,000 kilowatte, respect- 
ively, Рог Wr ОҒ 100,000 pounds and are proportionately increased. to 
approximately 500,000 pounds апа 100,000 kilowatts, а for Wy 
of. 150,000 pounds. 


Supersonic flight at а 30,000-foot altitude increases maximum 
reactor fuel-element temperatures to the 1200? to 1300? F range for approxi- 
mately the same values of airplane gross weight as required for subsonic 
flight; reactor heat releases required, however, are about four times 
greater than for subsonic flight. 


Supersonic flight at 50,000 feet altitude TEON nEs maximum reactor 
fuel-element temperatures to increase-to the 15009 to 15009 F range for 
Wy  0f.100,000 pounds, and to the 14000 to 1600? Е range for Мк of 
150,000 pounds for а reactor core diameter of-2 feet. For this flight 
Sondi tion, reduction in maximum fuel-element temperature to the 11009 to 
12009 F level for Мк of 100,000 pounds, and to the 12009 to 15009 Е level 
for. W. of 150,000 pounds, Gun Da achieved by ап increase in reactor core 
diameter to 2.5 feet. Airplane gross weight and reactor heat release are 
then 300,000 pounds апа 350,000 kilowatts, respectively, for Мк of 
100,000 pounds and are proportionally larger for Wr of 150,000 pounds. 


Enriched uranium investments for the hot unpoisoned reflected reactor, 
containing sufficient high-nickel-alloy fuel elements to provide appropri- 
ate heat-transfer surface for the aforementioned reactor and airplane 
flight conditions, are of the order of 35 and. 50 pounds for core diameters 
of 2.0 and 2.5 feet, resepctively. Ав fuel-element structural material in 
the reactor is reduced to его concentration, the uranium investments 
approach 15 and 20 pounds хок the 2.0- and 2.5-foot core diameters, 
respectively. The excess uranium required to counteract burnup and the 
poisoning effects of equilibrium xenon, samarium, and other fission 
products resulting from reactor operation for 24 hours at 300,000 kilowatts 
is estimated to be less than 10 pounds. At this poisoned condition, the 
temperature coefficient of reactivity for reactors with high-nickel-allay 
пет E TENERLE was calculated to be negative and of the order of -0.00006 
рег “Ж. 


INTRODUCTION 


Reactors cooled and moderated by liquid hydroxides, with uranium 
elther contained in fixed structural elements or present in compound 
form as а slurry in the hydroxide, have many attractive features to 
warrant detailed investigation of their applicability to aircraft nuclear 
propulsion. The high effectiveness of the hydroxides in slowing down 
neutrons makes for relatively small reactor sizes and hence small shield 
weights, which is particularly essential in aircraft application. The use 
of hydroxides, in functioning as combined coolant-moderator, leads to less 
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complicated reactor core structures compared with the air-cooled or 
liquid-met&l-cooled reactors which require incorporation of a separate 
moderator into the reactor. No pressurization is required to keep the 
hydroxides in the liquid state at the operating temperatures required 
for the aircraft applications; in addition, the hydroxides maintain 
reasonably high densities at these temperatures. These two character- 
istics represent important advantages over the use of water as a 
coolant-moderator. 


Important disadvantages in the use of the hydroxides are associated 
with: (a) their relatively low heat-transfer coefficients as compared 
with the liquid metals; (b) their high melting points; (c) their suscept- 
ibility to radiation decomposition; (d) their corrosive action on the 
available materials having the high-temperature strength required in the 

reactor and reactor coolant loop. 


The inferior heat-transfer ability of the hydroxides compared with 
the liquid metals is partly compensated by the greater reactor coolant- 
flow area and fuel-element heat-transfer surface area attainable with the 
hydroxides due to absence of a fixed moderator. Lower heat fluxes can 
therefore be attained in the hydroxide reactors so that temperature 
differentials between fuel element and coolant for the hydroxide case 
can be made to approach that for the liquid-metal-cooled reactor. 


Because of the high melting points of the hydroxides, the procedure 
for filling and for draining the reactor coolant requires prior heating 
of the reactor and loop, which introduces complications in starting up 
and shutting down the nuclear power plant. 


Whether radiation decomposition of the hydroxides at the operating 
conditions is an important enough consideration to rule out the hydroxides 
as possible reactor coolant-moderator remains unanswered at this time. 
Experience with water at high temperatures and pressures, has indicated 
that decomposition is much less severe than previously believed because 
of the high rates of recombination. Because the same mechanisms for 
recombination are operative for the hydroxides, this fact is encouraging. 


The major problem in the use of the hydroxides is the development of 
materials possessing both high-temperature strength and corrosion resis- 
tance to the hydroxides. Inasmuch as strength and corrosion-resistant 
properties of materials become progressively poorer with increases in 
temperature, it is advantageous to operate at the lowest Puel-element 
temperatures consistent with the requirements of reasonable reactor total 
heat releases and airplane gross weights. 


It is the purpose of this report to establish whether the hydroxide 
reactor shows sufficient promise for aircraft application from the stand- 
point of: (1) achieving small reactor sizes with reasonable fissionable 
material investments; and (2) providing required reactor heat releases at 
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required coolant temperatures without the necessity of excessive fuel- 
element temperatures. The investigation. was conducted. at the NACA Lewis 
laboratory. | 


Some evaluation work of the sodium hydroxide reactor for aircraft. 
propulsion Нав been made іп references 1 to 5. ОҒ the common hydroxides, 
sodium hydroxide has been generally selected because of its relatively 
low melting point and its satisfactory neutron absorption properties. 


Reference 1 is concerned with a low-powered stationary reactor 
intended to serve ав а step in the approach. to the high-power reactor. 
In reference 2 a general design and performance study is made of a 
nuclear-powered subsonic airplane having a homogeneous sodium hydroxide 
reactor (reactor in which the uranium is assumed as being dissolved or 
present as а slurry in the hydroxide). The study points up many of the 
practical problems involved in & hydroxide reactor power -plant system. 
In reference 5 а brief study is made-of a sodium hydroxide reactor with 
fixed fuel elements. The results indicste that the sodium hydroxide 
reactor with fixed fuel elements is sufficiently promising for subsonic 
aircraft propulsion (Mach number, 0.8) to justify developmental work. 


In the present report more extensive exploratory calculations are 
made of the criticality and heat-transfer characteristics of the sodium- 
hydroxide-cooled and moderated reactor with fixed fuel elements. The Р 
calculation results provide: 


(а) A generalized chart permitting rapid engineering evaluation 
of criticality requirements for reflected sodium hydroxide reactors 
incorporating structural materials of а wide range of M tL and 
concentrations. 


(b) An evaluation of maximum reactor fuel-element temperatures, for 
а range of reactor heat-transfer variables, necessary to satisfy the power 
requirements corresponding to both subsonic and ЕЕ aircraft 
propulsion. tet, a tees a TET Oo m u 

(c) A basis for compromise of the advantages of higher turbojet 
cycle efficiencies resulting from operation at high fuel-element tempere- 
tures with the advantages of higher structural-material strength and 
corrosion resistance resulting ый at low fuel-element 
temperatures. EE 1 g 


(а) An indication of the temperature coefficients of reactivity and 
excess uranium requirements corresponding to equilibrium fission-product 
poisoning for several reactor assemblies of interest. 
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SCOPE OF INVESTIGATION OF REACTIVITY VARIABLES 


The reactor core is taken ав а right circular cylinder of length- 
diameter ratio equal to unity. Fissionable material is contained within 
tubes or plates over which the sodium hydroxide NaOH flows to pick, up the 
fission heat generated within the fissionable material. The core is 
reflected by routing of МаОН around the core prior to its passing through 
the core, ав schematically indicated in figure 1. The composition of the 
reactor core is primarily NaOH with a relatively. small volume of structural 
material containing fissionsble material and possibly a diluent, or 
carrier, for the fissionable material. 


Criticality Calculations 


Reactor size апа uranium investment are dependent, to а large extent, 
on the amount and composition of the structural material present in the 
reactor core. At present, however, no structural material is known that 
is satisfactorily resistant to corrosion by NaOH at the temperatures 
required for the aircraft application. Hence, no definite assignment of 
structural material composition can be made for а specific evaluation of 
the NaOH-cooled and. moderated. reactor. 


As indicated in the section under "Method of Generalization of 
Criticality Results", reactor criticality requirements are influenced 
by structural material content principally through two quantities: the 
macroscopic thermal absorption cross section of the structure 2°, th 
and the reactor volume fraction of structure Ра. For опто ерее, the 
macroscopic cross section 28 A,th is put on & unit reactor volume basis 
and is hereafter referred to аб the thermal absorption parameter 
f A, the (Symbols are defined in appendix A.) The effect of variation 
of fo on criticality is determined on the basis of theoretical consider- 
ations. The effect of ғ” th for a given. f, must be determined by 
detailed criticality calculations made for several structural-materlal 
compositions and concentrations. Ву means of these parameters, criticality 
results for specific structural material contents in the reactor can be 
plotted in a generalized manner to be applicable for а wide variety of 
other structural material compositions and concentrations. 


The following specific reactor-core compositions were chosen for 
the criticality calculations to provide the basic data required for mapping 
the criticality characteristics of the family of NaOH-cooled and moderated 
reactors over а suitably large range of values of the thermal absorption | 


parameter £2", th. 
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0.0154. 
0.0074 


For these reactors, а wide i of uranium U concentration, 
corresponding to ratios of H to atoms R from 50 to 500, were 
investigated. Тһе reactor fission spectrum ranged from thermal to inter- 
mediate (fast-fission contributions from 10 to 50 percent of the total). 
Various thicknesses of NaOH reflector from zero to effectively infinite 
were investigsted. Reflector effects were measured by tbe two-group 
reflector savings, defined &s the difference ағына the unreflected and 
the reflected reactor-core radii. 


ln reactors I апа II, structural material content is taken ав 
8 percent by volume of Ni and Fe, respectively. Nuclear properties of 
these metals bracket the properties of a large number of high-temperature 
alloys апа ceramics considered to be satisfactory, nuclearwise, for use 
іп aircraft reactors. The 8-percent concentration of these metals is 
sufficiently high to provide for & fuel-element deslgn and auxiliary 
structure having adequate heat-transfer surface Рог а reactor of the 
order of 300,000-kilowatt output. Metallic Ма is included in the 
reactor as representative, nuclearwise, of the fissionable-material 
diluents or carriers that may be present within the fuel elements. 
Reactor ТЕГ represents the limiting case of zero structural material 
content. Reactor IV provides а comparison of the effect of the Na 
additive on reactor criticality and represents the truly homogeneous 
reactor type. 


Two-group neutron diffusion theory was employed in the criticality 
calculations. Spherical geometry was assumed in the reflected reactor 
calculations and & transformation made to cylindrical geometry under the 


assumption that the reflector savings remain the same. The derivation and 


method of solution of the two-group equations are briefly reviewed in 
appendix B. The evaluation of the gross nuclear constants required in the 


use of the two-group theory equations is briefly described in the following 


section "Evaluation of Nuclear Constants" and is presented in detail in 
appendix C. 


Results of the criticality calculations are presented in the 
following manner: 
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(1) For reactors I to IV, curves of the enriched uranium invest- 
ment as & function of reactor core diameter are presented for various 
reflector thicknesses. 


(2) Reflector savings for these reactors аге compared. 


(5) Representative neutron flux апа power generation distributions 
for reactors I and III аге presented. 


(4) For the generalized family of NaOH reactors, a chart of 
uranium investment ав а function of the thermal absorption parameter 


fel, th is presented for various reactor diameters and values of fg. 


Evaluation of Nuclear Constants 


Validity of the two-group representation used here depends 
entirely on the correctness of the nuclear constants used to represent, 
in a bulk fashion, the various: competing nuclear processes occurring 
within the two neutron-energy groups. Іп order to reduce the uncertain- 
ties involved in the evaluation of these bulk nuclear constants, the 
procedure used is patterned, as closely as possible, after that success- 
fully used in reference 4 to predict the criticality of small hydrogen- 
moderated (water, in this case) thermal reactors. The procedure for 
evaluating the constants, which is described in detail in appendix С, 
is briefly outlined as follows: 


А qn, f» 2д,?? Бр,» Pith: - The’ constants А qm, ғ Ё Af? Lf and 
Ph аге obtained by weighting local values according to the energy 
distribution of neutron flux in әп infinite medium of the same composition, 
as indicated by age theory. The fission neutron-energy spectrum (taken 
from reference 5) is included in determining the energy distribution. 


т^». - The solution of the transport equation for the second moment 


of the spatial distribution of neutrons siowing down from & point source 
іп an infinite medium consisting of а mixture of hydrogen and heavy 
elements, ав derived in reference 6, is used to determine ТРр. From 

the solution of reference 6, 2р із obtained as a function of fission 
energy and is then weighted over the fission neutron-energy spectrum. 

The calculations are normalized (see appendix C) to the experimentally 
determined value of 12», for water at room temperature, of 55 square 
centimeters. Тһе value of Lp thus obtained assumes no neutron : 
absorption during slowing down; for the reactors under consideration, 
however, significant absorption occurs in the energy range from about 
1000 ev to thermal. Inasmuch as the slowing-down process in this energy 
range can be considered continuous, age theory is approximately applicable 
and hence is used to estimate the effect of absorption on slowing down in 
this energy range. The difference in age with and without absorption in 


J- 


8 таса NACA НМ Е52119 


the energy range 1000 ev to thermal was calculated and found to be 
negligible. On this basis and, in addition, becatise the bulk of the 

Іл} contribution is due to the slowing down in the energy range above 
lOOO ev, the effect of$ absorption on Lp is taken to be negligible for 
the reactors under consideration. (Тһе effect of absorption in the fast 
group on criticality is hence accounted for solely through its effect on 


Рен). 
ATR, һе - The method of reference 7 for estimating the effect of 


2525 


chemical binding on thermal neutron diffusion.in а hydrogenous medium 

18 used to evaluste ATR, th- In this method the experimentally determined 
variation of scattering cross section og, with neutron energy for hydrogen 
(ав measured for water, reference 8) is used to obtain Жв ав а function 
of neutron energy. The assumption involved here is that the chemical bond 
of H in NaOH is the same as that in water. The value of Атв,ър 16 
then obtained by weighting App according to the flux of neutrons in a 
Maxwellian distribution corresponding to the moderator temperature. 


WEE Ly th. - Local values of Za and Бр are assumed to vary 


with neutron energy according to 1/у in the vicinity of thermal energy; 
5 th and 5а th аге then obtained by weighting Z4 and Фр, respect- 
lvely, &ccording to the fiux of neutrons in & Maxwellian distribution 
corresponding to the moderator temperature. 


The constants Кр, Кер, and Loan are then given by: 


viz fr > Е th ^ 

р 2 2 TR,th 
Ko = кое H К+ = po 2 ds = — 
5 A,f у n A th Su 3 А, th 


In the criticality calculations, the average temperature of the 
NaOH moderator is taken equal to 14500 F (Еру = 0.092 ev). In all the 
reactors, the constituents of the reactor core are assumed to be completely 
intermixed so that the criticality calculation results до по include any 
self-shielding effects arising from any heterogeneity of the reactor 
composition. The fuel cross-section data were taken from measurements on 
K-25 end product consisting of 91.5 percent 0855, 1.5 percent 0254, ара 
7.0 percent 0258 normalized per atom of 0250. Densities of the materials 
in the reactor (at 14509 F average temperature) are tabulated in table I. 


A summary of two-group constants for the specific reactors analyzed 
in the present study is presented in table II. 
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Calculations of Temperature Coefficient of Reactivity 
апа Excess Uranium Requirement 


In order to provide an indication of the static stability 
characteristics of NaOH-cooled and moderated reactors, two-group per- 
turbation theory was used to calculate the temperature coefficient of 
reactivity for the following cases: 


Reactor I: R= 100 and 200, ty = 6 inches 
Reactor III: R= 100 and 400, t, = 6 inches 


First-order perturbation formulas, modified to include the effect 
of fast fission, are presented in appendix D. 


Calculations were made for the hot reactor with equilibrium Хе155 
and 24-hour 51289 апа other fission-produet poison concentrations 
corresponding to reactor operation at 500,000 kilowatts; these concen- 
trations were determined on the basis of information presented in refer- 
ences. 9 апа 10. The poisons are assumed to be distributed uniformly. 
over the reactor-core volume. 


Contributions to the temperature coefficient of reactivity 
included herein arise from: (1) The change in neutron energy range 
constituting the fast group, (2) the change in density of Маон (and 
hence in atom density of Ма, О, апа Н), and (3) the change in thermal 
microscopic cross-sections. 


The linear variation of density of Маон with temperature was 
taken from reference 11. With the exception of Хе155 ana 5149, 
absorption and fission microscopic cross-sections in the vicinity of 
thermal energy are assumed to follow the 1/у law. The variation of 
microscopic absorption cross-section of Xel55 with temperature, as 
averaged over the Maxwellian distribution, is obtained from reference 9. 
The variation of microscopic thermal cross-section of 81149 with 
temperature (which contribution to temperature coefficient of reactivity 
is small relative to the Хе155 contribution) was estimated on the basis 
Of the даба of reference 10. 


Appendix D presents the perturbation formula and procedure for 
obtaining the adjoint functions used in the formula; the perturbation 
formula, presented is an extension of that derived in reference 12 to 
include the effects of fast fission. Appendix D also presents details 
of the method of evaluation of the various weighting factors contributing 
to the temperature coefficient of reactivity. 


In addition to determination of temperature coefficient of 
reactivity, two-group perturbation theory was used to estimate the 
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excess fissionable material required to compensate for burn-up and 
equilibrium fission-product poisoning for the same reactors. Although 
the principal poisons are Хе! 59 and Smi49, s host of other fission 
fragment poisons appéar and are accounted for by the assumption that for 
every uranium atom fissioned, the equivalent of a single atom appears 
having а microscopic absorption cross section of 100 barns at 0.025 ет. 
The details of the method of calculation are outlined in appendix D. 

The results of the aforementioned calculations are presented in tabular 


form. 


Method of Generalization of Criticality Results 


For the range of NaOH reactor compositions of interest for aircraft, 
the principal contribution to neutron slowing-down and diffusion processes 
is made by the NaOH, whereas the principal contribution to the absorption 
process is made by the structural material and uranium. For engineering 
evaluations, therefore, it may be assumed that the over-all scattering 
and slowing-down in the reactor core are relatively unaffected by 
variation in concentration. and composition of structural material from 
the values assumed. for specific reactors I to IV. 


Therefore, the only two important effects on reactor criticality 
requirements due to variation in structural material composition and 
concentration in the reactor аге: (a) the effect of neutron absorption 
by the structural material, which is given by the thermal absorption 
parameter 70, 4, for а l/v absorber, and (b) the effect of dis- 
placement of NaOH moderator by the structural material when volume con- 
centrations fg differ from the values assumed for the specific 
reactors. pe з мен 7 mE | 


Variation in Ға about а reference value can be treated as 
equivalent to the effect of void space; а decrease in structural-material 
volume from the reference value is equivalent to the removal of void 
Space by the same amount; an increase іп structural-material volume from 
the reference value is equivalent to displacement of moderator by the 
same amount. Тһе effect of void врасе on reactor size and uranium invest- 
ment, that is, the effect of fg, is treated in appendix B. 


Ás previously mentioned, the effect of ғ.” А + 18 obtained from 
specific criticality calculations in‘which fgi5, 4, is varied over the 
range of interest. For a mixture of elements, the parameter fg254 th 


ів given by: 
В. 5 ~ ii | M 120% | 
Pod A,th.= leg 2л,1 T Со Za ,2 г... + CH Zain) (1) 


4 - - PL = 


where c- ів the fraction by volume of elements in structure. 
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In table IIT are listed representative alloying elements which may 
be employed. These absorbers follow the 1/у law in the thermal region 
with little probability of strongly absorptive, low-lying resonances 
being present (see reference 13). Therefore the generalization may also be 
used to evaluate criticality requirements for NaOH reactors containing 
materials other than Ni and Fe and for a range of volume concentration. 


Aithough the generalization is not strictly rigorous in that it 
neglects the small differences in neutron slowing-down properties due 
to other structural materials, it does account for the important first- 
order effects provided: 


(1) The neutron absorption cross section of the structural material 
varies with neutron velocity according to 1/у at least up to energies 
of the order of 100 electron volts. 


(2) The volume of structural material in the reactor is small 
compared with the volume of moderator. 


SCOPE OF INVESTIGATION OF PERFORMANCE OF NUCLEAR-POWERED ATRCRAFT 


Ihe purpose of this performance study is to determine the airplane 
gross weight and reactor heat-release requirements as affected by the 
maximum temperature of the reactor fuel elements for & range of reactor 
and reactor heat-transfer variables for various airplane flight conditions. 


Ihe turbojet cycle, involving а tertiary system of two separate 
closed liquid circuits and ап open &ir cycle, is used in the study of the 
performance of а nuclear-powered aircraft utilizing the NaOH reactor as 
the heat source. А schematic diagram of this cycle is shown in figure 2. 
In the primary liquid circuit, NaOH is heated as it flows through the 
reactor and then cooled as it flows through the primary heat exchanger 
where it gives up its heat to а liquid metal flowing in the secondary 
liquid circuit. Іп the secondary circuit, heat picked up by the liquid 
metal is transferred to the air in the secondary hest exchanger. 


Air enters the diffuser of the turbojet engine, is compressed by 
а compressor, is heated in the secondary heat exchanger, and is expanded 
through а turbine which extracts sufficient energy from the high-pressure, 
high-temperature air to run the compressor. Finally the air expands 
through an exhaust nozzle into the atmosphere to provide the propulsive 
jet thrust. 


The configuration of reactor fuel elements chosen Рог determi- 
nation of heat-transfer characteristics of the NaOH reactor is shown 
schematically in figure 1. This configuration was chosen for convenience 
in the heat-transfer study intended to indicate what may be accomplished 
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with respect to beat transfer іп the NaOH reactor; no detailed study of 
the structural, mechanical, or fabrication problems associated with 
various core designs was made in arriving аб the configuration chosen. 


As indicated in figure 1, the fuel elements are considered to be 
plates of.sandwich-type construction arranged in a closely spaced parallel 
array within the reactor core. NaOH flows in a single pass between the 
plates from one face of the cylindrical core through to the opposite face. 
(Choice of small-diameter tubes instead of plates would give approximately 
the same heat-transfer results.) 


2529 


Reactor Heat Transfer 


In the-heat-transfer study, the maximum temperature of the fuel- 
element plates corresponding to any value of total reactor heat release 
and reactor-exit coolant temperature is related to the variables: reactor 
core diameter D a(note, Le = Пе)» fraction of core volume occupied by 
fuel-element шас ier fg» fuel-element plate thickness Up» and coolant 
velocity У. 


Two methods of reactor operation have been assumed: (1) uniform 
heat generation over entire reactor-core volume, and (2) uniform fuel- 
element wall temperature. 


Basic relations and procedure used in the heat-transfer study— 
are outlined as follows: 


(1) From geometrical considerations, the plate spacing sp, area 
for coolant flow А, and heat-transfer surface area 5 are obtained for. 
assigned reactor core size Па, percent of соке volume occupied by fuel- 
element material fg, and plate thickness өр. 


(2) Тһе heat-transfer coefficient h corresponding to a given 
coolant velocity V is then obtained by use of the Nusselt relation 


0.8 0.4 
h De (ша сри 
k : " k 


The physical properties for the NaOH coolant employed in the 
present calculations. аге as follows: 


where De ж 2 Бр. 
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Density, ps 1b/eu ft е е * • . е Ф ө е . Ld а . . * . . е a е ~ 98 „© 
Specific heat, ст, Btu/lb-9F . E сын б тв а doped asse БС. 
Viscosity, І, lb hr-ft • • • . . в е • е е ә е • . s . 2.5 


Conductivity, k, Btu/(hr)(sq ғы)(оғ/ға). Мир и Ре Ко 


Density data were taken from reference 11. Viscosity, specific heat, and 
thermal conductivity of NaOH were obtained from даља transmitted from 


Battelle Memorial Institute. 


(3) For the case of constant rate of heat input along а flow 
passage, the temperature differential between the passage wall and coolant 
is constant along the flow passage (for constant h along passage) я 
Hence the heat balance пау be expressed as: 


H = bs Са 7 Te ex) = pAVe,, (To ех ә Те еп) "S E 


TO -T T -T P 
е ее апа == аге readily M 
determined. Тре temperature differences аге expressed per kilowstt of 
heat release in order to make the heat-transfer results general for any 


value of reactor heat release. 


from which the ratios 


(4) For the case of uniform wall temperature, the temperature 
differential between the passage wall and coolant varies along the flow 
passage. The heat balance for this case is written as: 


(те P ех-Те P en) 


H = hS = РАУср (Те, ex-Tc , en) 


О 
T Ws ex te, en 


loge 
Ws ex Іс „ех 
О 
Т w,ex Te, av Te,ex"tc,av 
from which the ratios е: зава and --- are determined. 


(5) The temperature drop across the fuel-element plate is calcu- 
Ws ex 2, ex 

2 2 
assuming all the heat is generated at the center line of the plate cross 
section. The following relation can then be obtained: 


lated by use of the heat-conduction equation u/s м - Ка 


Tw, ex tc, av " Ти, ex-Te av + tp 
H 7 H 2k.S 
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The thermal conductivity of the fuel plates used in calculation of 
temperature drop across the plates is teken as 28 Btu/(hr)(sq ft)(?rF/ft) 
the value for nickel at 1500 Е. The thermal conductivity for carbon, 
stainless, and high-alloy steels is about 15 Btu/(hr)(sq ft)( F/ft). 
(reference 15).  Inasmuch ав the temperature drop across the fuel plates, 
for the conditions considered, is generally less than IO percent of the 
over-all drop to the average coolant temperature, variation in thermal 
conductivity of the fuel-plate material has not been considered. 


The results of the foregoing procedure аге presented in the follow- 


ing manner: "v т E 


For & fixed reactor core size and percent of core volume occupied 


T -т T -T 
by fuel-element material, eee and = are plotted 


against V (range, 5 to 50 ft/sec) for various fuel plate thicknesses. 


Reactor heat-transfer characteristics are presented for the case of 
uniform heat generation for: 


Reactor core diameter, ft. e ..:...:.... e я o 2.0, 2.5, 3.0 
Fuel-element reactor volume fraction . . .. ........ 0.06, 0.12 
Fuel-element plate thickness, in. . . . e . è . > . 0.012, 0.016, 9.020 


Comparative heat-transfer characteristics are presented for the 
case of uniform wall temperature for each of the reactors Рог а fuel- 
element-plate thickness of 0.012 inch. | | 


Airplane апа Turbojet Cycle 


The performance study of the turbojet cycle system in а nuclear- 
powered aircraft 18 based on the optimized engine performance results .. 
of reference 14, wherein the turbojet cycle 1в optimized to give minimum 
gross airplane weight Wg for fixed values of airplane lift-drag ratio 
L/D, structure-to-gross-weight ratio "Не, апа weight of shield, 
reactor, pay load and auxiliary equipment Wr. The optimized performance 
is given by relations between the air heat-exchanger effective wall 
temperature T, app and values of the following variables: (а) engine 
thrust per weight of engine plus air heat exchanger Рат, (b) net 
thrust per pound of air per second Fn/Was and (c) heat addition per 
pound air Ah,. From these quantities, the airplane gross weight, 
reactor heat release, and engine air flow are found from the following 
relations for any value of air heat-exchanger effective well temperature: 


Wy | 
е Ws + 


Ча (тым) (1/0) 
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Н = Wa Ah, 


Values of variables selected for the airplane study are: Ук, 100,000 
and 150,000 pounds; L/D = 6.5 (supersonic) or 18 (subsonic); and 
Ws e = 0.50. 


The values ОҒ Wr selected are believed to be representative of 
the divided and unit-type shields contemplated for the nuclear airplane. 
The values of L/D selected conform with existing practice for the sub- 
sonic case and represent a reasonable value for the supersonic case. 

The value of Ws/We used conforms with the findings of recent nuclear 
airplane weight analyses. | 


Maximum Fuel-Element Temperature Evaluation 


in order to relate the airplane gross weight and reactor heat 
release to the maximum temperature of the reactor fuel elements, it is 
necessary to combine the engine performance and the heat-transfer studies. 
The procedure used is outlined as follows: 


(а) For assigned values of Wx, L/D, апа W,/W,, the method out- 
lined in the previous section is used to determine 02, Was and Н for 
а range of values of Ту, еРР · 


(b) On the basis of а preliminary primary heat-exchanger design 
study, the average NaOH coolant temperature within the reactor is taken 
as 100° F higher than Tw, eff: that is, To ay - Te, err = 1009 Е. The 
heat-exchanger design study indicated that this temperature differential 
can be easily attained with reasonably small primary heat-exchanger size 
for reactor heat releases required for supersonic flight conditions. 


(c) For assigned values of reactor size, coolant velocity, fuel- 
element material volume, and fuel-element thickness, the maximum fuel- 
element temperature is computed by adding the following temperature 
differences: 


Ту,ех = (Т ех E То ау) * (To, av 3 To erf) 


А schematic diagram of the heat-exchanger system and relative 
temperature distribution is presented in figure 3. The temperature 
difference Ty,ex - Tc,ay 18 found from the heat-transfer parameter 


Ty 9 ex" Te ‚ах 


Я — by multipiying by the value of' H determined in step (а). 
есу 
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RESULTS AND DISCUSSION 
Specific Reactor Calculations 


Criticality calculation results for the core compositions 
represented by reactors I to IV are presented in figure 4. Іп these 
figures, uranium investment is plotted against cylinder core diameter 
for reflector thicknesses of 0, 5, 6, and 12 inches. (Cylinder length- 
diameter ratio is taken equal to unity.) Іп the calculation of these 
curves for each core composition and reflector thickness, а range of 
' values of atom ratio В is taken and core diameter and resultant uranium 
investment for criticality are determined. Lines of constant В are 
given on each figure. 


2585 


The curves for each reactor-core composition indicate а minimum 
total uranium investment for each reflector thickness; minimum total 
investment decreases and occurs at progressively smaller core diameters 
as reflector thickness increases to effectively infinite values. 


The existence of а reactor size for minimum total uranium investment 
may be explained &s follows: Ав uranium concentration (per unit volume) 
is increased, the neutron-leakage tolerable for criticality is increased; 
а decrease in core diameter results. 


At the large соке diameters, the-neutron leakage is relatively 
lnsensitive to diameter so that in order to maintain criticality, а 
large percentage decrease in core diameter occurs for a small percent- 
age increase in uranium concentration; hence, the total uranium invest- 
ment, which is proportional to the uranium investment and the cube of 
the diameter, decreases. At the small core diameters, the neutron leak- 
age is extremely sensitive to diameter and so the opposite effects occur. 
At some intermediate value of core diameter, the decrease in diameter 
exactly counteracts the increase in uranium concentration with respect to 
total uranium investment; this defines the core diameter for minimum 
investment. 


Minimum total uranium investments and corresponding reactor-core 
diameters for unreflected and reflected reactors are summarized from 
figure 4 in the following table: | 
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Reactor Composition by Minimum uranium | Reactor core diameter 
yolume investment - (ft) 
(1b) 
NaOH | Ni Unre- 6-in. Unre- 6-in. 
flected|reflector | flected | reflector 

I 0.20 95 52 2.7 1.7 
II .82 69 25 2.7 du 
III .90 43 15 2.7 3..5 
IV 1.00 58 11 Рио 1.5 


Considerations of reactor heat transfer and high-altitude supersonic 
Plight power requirements, presented later, indicate that core diameters 
of the order of 2.5 feet are required to avoid unreasonably high fuel- 
element temperatures. For a 2.5-foot core diameter and a 6-inch 
reflector, figure 4 indicates investments of 51, 36, and 21 pounds for 
reactors I, II, and ТТТ, respectively. 


Reflector Savings 


The variation of NaOH reflector savings (defined as the difference 
between the unreflected and reflected reactor core radii) with reflector 
thickness for reactors I to IV is presented in figure 5 for representative 
values of atom ratio R. Reflector savings increase markedly for reflector 
thicknesses up to 6 inches, but level off rapidly for larger thicknesses. 
А line indicating reflector savings equal to reflector thickness is 
included in figure 5. Reflector savings are about equal to reflector 
thickness up to thicknesses of 3 inches, indicating no net change in 
reactor-core radius-plus-reflector thickness; for thicker reflectors, 
reflector savings are smaller than reflector thickness so that the over- 
all core-plus-reflector dimensions become larger. From considerations 
both of over-all size (core plus reflector) and of reflector effective- 
ness in reducing uranium investment, а reflector thickness of about 
6 inches appears to be а satisfactory compromise. 


Ine different reactor core compositions considered herein exhibit 
approximately the ваше reflector savings for any given value of reflector 
thickness. The small effect of increased nonproductive absorber in the 
core on the reflector savings may be noted in the slightly reduced 
reflector savings attained for reactors I and II compared with reactors 
ТТТ апа IV. Variation of В іп the core also has little effect on the 
magnitude of reflector savings. 
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Although the magnitudes of the savings are about the same for 
thermal reactors (Тож concentration of uranium) for which соге diameters 
are large, апа for intermediate reactors (high concentration of uranium) 
for which core diameters are relatively small, the savings аге mueh 
smaller percentages of the core radii for the thermal reactors because 
the neutron leakage is smaller for thermal reactors. 


The curves of figure 4 for the 6-inch NaOH reflector thickness are 
replotted in figure.6 for convenience of comparison of reactors I to IV. 
In every case, the reactor core diameters corresponding to minimum invest- 
ment sre less than 2 feet. Тһе curves become steeper on both sides of 
the minimum investment point as the nonproductive absorption is progress- 
ively increased from reactor ТТТ to II to-1; this steepness is more pro- 
nounced in the region of core sizes below that for minimum investment. 
Inasmuch as the uncertainties introduced by two-group theory become 
smaller as core size is increased (the value of R increased with the 
result that fast-fission contribution was decreased), greater confidence 
can be placed in the results for core sizes to the а а of the minimum 
investment values. - - 5 = 


Criticality Generalization 


Chart construction. - On the basis of the considerations presented 


in the section entitled "Methods of Generalization of Criticality Results", 


the specific results for reactors I, II, and ТТТ are replotted to be 
applicable, in engineering evaluations, to & wide variety of structural- 
material concentrations and compositions, subject to the restrictions: 
(a) that the material does not exhibit any large neutron cross-section 
resonance near thermal energy (below about 100 ev), and (b) that the 
concentrations must be relatively smell (say, less than 16 percent). 

The generalization, results of which.are plotted in figure 7, is 
accomplished as follows: 


А plot (fig. 7(&)) is made of uranium investment М9 against the 
thermal absorption parameter РБА, th Ber various unreflected розе 


core diameters conposed ав follows: 


Volume fraction of strueture,. e e . e s e è e © e e e 2а 
Volume fraction of NaOH plus Na, 5 в . е є е е s e e в. е e в. е . 1 - fg 
N&/NaOH volume ratio,. іт WO ка = = =. ab а за. с 122 


The Specific results for reactors I, ІІ, апд III аге used in the 
construction of this plot employing а reference value of Ға of 0.08. 
For reactors I and II, the value of Ға is 0.08 and so the results are 
directly applicable; for reactor III, fg is zero so that & void correction 
to fg of O.08 was made by the method discussed in appendix B. 
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The unreflected core ш indicated on figure 7(а) аге referred 
to ав the chart diameters (= 35 По. The slant lines on the left of 
this plot yield the unreflected reactor total uranium investments МУ 
for any assigned value of fg. 


The insert plot (fig. 7(b)) indicates the reduction in core dia- 
meter possible with various thicknesses of NaOH reflector. The 
reduction in core diameter is equal to twice the reflector savings as 
presented in figure 5; inasmuch ав reflector savings are relatively 
insensitive to core composition, average values were selected for 


figure 7(Ъ). 


Therefore, for assigned values of fg and По of interest, 
figure 7(а) enables determination of enriched uranium investments for 
the unreflected reactors for any value of Ё} л ақ. Figure 7(b) permits 
the evaluation of the reflected core diameter for any v&lue of reflector 
thickness; the uranium investment for reflected reactors, which is pro- 
portional to the cube of the core diameter, may be calculated from ihe 
uranium investment for the unreflected core. 


Chart procedure. - The procedure for using the chart is given by 
the following steps: | 


1. A reflected NaOH reactor with.actual structural material con- 
centration fg and core diameter D,(t,) is assigned. 


2. Figure 7(b) gives the reduction in reactor diameter due to NaOH 
reflectors of various thickness; tbe unreflected core diameter of interest 
is then evaluated. ` 


De = Deltr})} + AD,(t,) 


-f 
3. The unreflected chart reactor diameter Е = 


5 33 De is calculated 
by means of the actual value of fg and the unreflected core diameter of 


interest Do- 


4. Figure 7(a) yields а value of uranium investment М9 for the 
unreflected reactor with the actual value of Ға. 


5. The uranium investment for the reflected reactor of interest 
is then computed. 
| 5 
ud 


W(t.) = W ск 


The use of the chart can best be illustrated by an example: 
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Example. - It is desired to find the enriched uranium investment 
for а cylindrical reactor of-length-diameter ratio of unity, moderated 
by NaOH (with small percentage of Ма additive) апа reflected by а 
6-inch thickness of NaOH. The reactor core is 2.5 feet in diameter and 
is to contain volume concentrations of Inconel of (a) 8 percent, and 
(b) 12 percent. The thermal absorption parameter foD th 18 evaluated 


for each case іп the following table by equation (1): 


Element| Fraction |Fraction M Case & | Case b 
by weight|by volume|(from Е АВ nn i Шр, ВОД, ЪН е РА th 


0.12 X 
0.159 


0.08 X 


Сто | оо | = | $359 9:38. 


The procedure for using the general criticality results of figure 7 
to evaluate the unreflected and reflected reactor core diameter and 
uranium investments 18 indicated in the following table: 


р.(%,.) Аре AA 


0.08) 0.0127 | 2.50 | 0.76| 5.26 
. 12 .0191 2.50 .76| 5.26 


The generalized results apply for reactors operating in the thermal 
temperature range of.14009 to 15009 Е; for lower temperatures, investments 
are Slightly lower than indicated by figure 7. 


Typical Neutron Flux and Heat-Generation Distributions 


Two-group neutron flux ара heat-generation distributions have been 
determined for reactors I and III and аге presented herein ав representa- 
tive of NaOH-cooled and moderated reactors. The distributions are pre- 
sented as а function of cylinder radius for an equivalent spherical 
reactor as an indication of the distributions for the cylindrical 


geometry. 


Neutron flux distributions. - The fast and thermal neutron flux 
distributions are shown in figure 8 for the 6-inch reflector thickness 
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and for а value of atom ratio В of 100 (indicated in figure 4 to be 
near the minimum investment point). For convenience, the flux values 
shown correspond to reactor total power output of 500,000 kilowatts; 
inasmuch as fluxes are directly proportional to power, ihey may be 
adjusted for power outputs other than 500,000 kilowatts. Flux levels 
for the same total power output are higher for reactor III than for 
reactor I principally because of the smaller core size possible with 
reactor ТТТ. 


Because of the smaller absorption in reactor III, thermal neutrons 
resulting from slowing down in the reflector can penetrate’ farther into 
the core of reactor III than reactor I. Hence, the effectiveness of the 
reflector toward flattening the thermal flux in the core is greater in 
reactor III than I; for example, the ratio of minimum to maximum thermal 
flux in the соке is 0.65 and 0.45 for reactors ТТТ and I, respectively. 


The ratio of fast to thermal neutron flux is about 9 for reactor I 
and about 7 for reactor IIT. These ratios will decrease as reactor size 
and corresponding investment are increased and less dependence upon fast 
fissions is required. 


Heat-generation distributions. - Radial heat-generation distributions 
are shown in figure 9 for various values of reflector thickness апа for 
the value of R equal to 100. Figure 9 is а plot of the ratio of local 
specific heat release to average specific heat release against reactor 
core radius expressed as а fraction of total core radius. 


The distributions demonstrate the degree to which reactors with 
NaOH reflectors шау be expected to approach uniform heat generation with 
uniform distribution of uranium. The effect of reflector thickness 
saturates rapidly with little change in distribution above thicknesses 
of 6 inches. 


Larger variations in local heat release with radius are to be 
expected of reactors containing appreciable concentrations of absorptive 
structure, inasmuch as these reactors depend in greater measure upon the 
fast flux to provide fissions than reactors with little absorptive struc- 
ture. Іп addition, as previously mentioned, neutrons which enter the 
core after being slowed down in the reflector do not penetrate as far 
into the core for reactor I as they do for reactor III, thus contributing 
to further nonuniformity in the power-generation distribution. 


These two points аге iliustrated in figure 10 in which the fraction 
of total local fissions produced thermally is presented as a function of 
core radius for reactors I and ТТТ for various reflector thicknesses. 

It may be seen that reactor I has about 6 percent fewer thermal fissions 
than reactor III throughout the core. Іп addition, the high concentration 
of thermal fissions at the core-reflector interface falls off more rapidly 
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for reactor I.than for reactor III. It is noted that these reactors can 
be classified as intermediate reactors in that about 40 percent of total 
fissions are produced in the fast neutron group. 


The heat-generation distribution over the reactor core volume can 
be varied by nonuniform distribution of fissionable material. Refer- 
ence 16, which presents Ғіввіопар1е-шафегіз1 distributions corresponding 
to uniform heat generation for & representative thermal reactor assembly 
with various reflector thicknesses, indicates that appreciable improve- 
ment in power-generation distribution can be achieved in a well reflected ' 
reactor by redistribution of the fissionable material with only a small 
additional investment of fissionable material (of the order of 10 to 
15 percent for the assemblies of reference 16). Although the results of 
reference 16 are for thermal reactors in which the fissionable material 
was redistributed to give uniform heat generation, the same general 
conclusions should apply for the intermediate reactors presently under 
consideration for which the fissionable material may be redistributed to 
attain uniform fuel-element temperature SBCOHEROUC the reactor core 
volume. 


Static Stability Characteristics and Excess Uranium Requirements 


The—temperature coefficient of reactivity апа excess uranium 
requirements have been.computed for reactors I and ITI. Two гергевепба- 
tive examples of each reactor have been considered (an intermediate 
reactor of small core size and а more thermal reactor of relatively large 
core size) in order to illustrate both the effect of structure and the 
effect of neutron leakage on reactor stability. Reactors reflected with 
& 6-inch thickness of NaOH have been considered for both cases. 


The radial distribution of adjoint functlons Фе ала Pth 


relative to the value of P tth at the reactor axis, for the smaller 
hot, unpoisoned reactors I and III with an atom ratio R of 100, are 
presented in figure 11. These adjoint functions, computed in accordance 
with the methods described in &ppendix D, аге used in conjunction with 
the actual neutron flux distributions to provide weighting factors which 
evaluate the relative importance of small local changes or perturbations 
of the reactor constants at any reactor geometrical position. 


For example, the effect on pile reactivity of & small change in 
thermal absorption cross section Едуин at а particular position of the 
reactor core is proportional to the product of the thermal flux and the 
thermal adjoint function at that position. The over-all effect on 
reactivity of small changes in ША, th Over the entire reactor core is 
evaluated by the integrated effect over the entire core volume as shown 
in appendix D. 
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Temperature coefficient of reactivity. - Temperature coefficients 
Of reactivity have been computed for both poisoned and unpoisoned 
reactors operating at 14509 Е and а power of 300,000 kilowatts. Excess 
uranium requirements have been estimated for these reactors on the - 
assumption that the reactors were predominantiy thermal. The poisoned 
condition considers the equilibrium xenon concentration and concentrations 
of 24-hour accumulation of samarium and fission products for operation at 
а constant power of 500,000 kilowatts; these concentratiéns are listed оп 
the following table for the reactors considered: 


Reactor н/у . Ny Ма, Neiss. prod. 
atom ratio, уреди ( atoms/cc 


core core) 


Reactors are reflected by 6-inch thickness of N&OH. 


individual temperature coefficients of reactivity associated with 
changes with temperature of the two-group parameters for each of the 
four reactors are listed in table IV. Derivatives of each parameter with 
temperature were obtained numerically from the calculated variation of 
the parameter with temperature and satisfactorily represent actual 
derivatives for the temperature range from about 1200° to 1700° F. 
Temperature coefficients of reactivity ар/ат were computed by the 
relations given in appendix D. 


In the calculation of ар/ат for reactor III which does not contain 
fuel-element structure and would therefore physically correspond to the 
homogeneous reactor, it has been assumed that the uranium concentrations 
remain fixed and thet only the cross section varies with temperature; 
results for reactors І and III given in the table are then directly 
comparable. Of course, the reduction in uranium concentration accompany- 
ing the expansion of NaOH in a homogeneous reactor would produce & sub- 
stantially more negative temperature coefficient of reactivity for 
reactor ТТТ. | 


From table ТУ it may be noted that the combined contribution to the 
temperature coefficient of hath (without poison), Кен» ZA, th? LAP? 
and Ке БАР, that is, of the fission and absorption processes, is 


relatively small. The effect of the poisons (given by difference in 
contributions of 24th with and without poison) is slightly positive. 


aes 
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The combined contributions of ATR, ths Атв,р» and У that is, of 


the scattering and slowing-down processes, give а relatively large 
negative temperature coefficient which is associated mainly with the 
decrease in density of moderator with increase in temperature. Ав. 

а result of this large negative contribution, the over-all temperature 
coefficients of the reactors considered are negative. The magnitude of- 
the negative temperature coefficient increases as the size of the 
reactor decreases. 5 
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Net values of the temperature coefficient of reactivity ар/ат 
vary from about -0.00006 for the small-sized reactors to about -0.00002 
for the large reactors. А temperature coefficient of reactivity of 
about -0.00006 per °F results in а change in reactivity of -0.6 percent 
for а 1009 Е temperature rise; this reactivity is of the same order of 
magnitude: as that associated with the delayed-fission neutrons and во 
indicates the negative coefficient to be significant from the point of 
view of steady-state reactor self-control. 

\ 

Excess uranium requirements. - The excess uranium requirement 
AWU Рог the poisoned critical reactor over the investment required for 
the unpoisoned critical reactor. WÙ is also listed in table IV. It may 
be seen that the excess uranium requirement necessary to maintain eriti- 
cality for these reactors after 24&-hour operation at а power of 500,000 
kilowatts increases ав the size of the reactor producing this power 
decreases. The excess uranium requirement appears to be less than 
10 pounds for ali the reactors herein considered. 


Reactor Heat-Transfer Characteristics 


It is recalled that the reactor heat-transfer characteristics 
presented are independent of any airplane considerations. The results 
are presented per kilowatt of total reactor heat release and are 
related to the airplane only when the power requirements for & particular 
flight condition are established. 


Figure 12 presents plots of the difference between maximum fuel- 
element temperature and average NaOH coolant temperature per kilowatt 


of total reactor heat release, че ме, against coolant velocity. 


Each plot 16 Рог а given core diameter (2, 2. 5, or 5 ft) and concen- 
tration of structural material (6 or 12 percent by volume) for fuel- 
element plate thicknesses of 0.012, 0.016, and 0.020 inch. Heat gener- 
ation is assumed to be uniform over the reactor core for the 0.016- 
and 0.020-inch plate thicknesses; for the 0.012-inch plate thickness, 
both cases of uniform heat generation and of uniform fuel-element wall 
temperature are presented. Also included in each plot is the difference 
between the exit and average NaOH temperature. The effects of pertinent 
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variables оп the difference between maximum fuel-element temperature 

amd. .exit coolant temperature аге directly obtainable from figure 12. 
The following tabulation of results, obtained from figure 12 for а 


reactor with total heat release of 500,000 kilowatts, indicates the 
effects on Ту ех - То,ау ОҒ representative values of the heat-transfer 


variables. 
Reactor Ty, ех" te Р ат 
oe изе) Ge.) operation (OF) 
590 


“Uniform heat generation 
Uniform wall temperature 


Airplane and Turbojet Cycle Characteristics 


Flight condition requirements. - Figures 15 апа 14 present plots 


of reactor heat release Н, апа &irplane gross weight Ме (expressed 
in ratio to shield, reactor, pay-load, and auxiliary weight Wy) against 
effective wall temperature of the air heat exchanger Тұ, eff for flight 
at 30,000- and 50,000-foot altitudes at Mach numbers of 0. 9 and 1.5. 

The optimized turbojet сусје conditions determined in reference 14 for 
each of these flight conditions were used m the construction of these 
plots. 


Figure 13 shows that the effect of flight speed on reactor heat 
release is much greater than the effect of altitude, chiefly because of 
the greatly reduced airplane lift-drag ratios (and hence increased power 
requirements) encountered at supersonic speeds. For Tw, eff about 
17009 F, heat releases required for supersonic flight are of the order 
of three times greater than for subsonic flight; for "Ty eer about 1100? F, 
the heat-release requirements are of the order of five times greater. For 
the range of values of Т, epe plotted in figure 15, heat releases Рог 
subsonic flight are fairly insensitive to variation in Ty,epr; for 
supersonic flight, heat releases are quite sensitive to Ty epp, Parti- 
cularly in the range below 1200? Е. 


- 
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Figure 14 shows that reduction of altitude from 50,000 feet to 
50,000 feet for_the 1.5 Mach number significantly reduces gross weight . 
requirements to the point where the gross weights &re comparable with 
subsonic v&lues. 


Maximum reactor fuel-element temperature. - The reactor heat transfer 
апа the turbojet cycle characteristics have been combined to establish the 
magnitude of maximum reactor temperatures required for the four flight 
conditions considered. 


Maximum fuel-element temperatures аге plotted in figures 15 and 16 
against reactor heat releases required for the four flight conditions 
considered, for reactor operation st uniform heat generation апа uniform 
wall temperature, апа for the following conditions: 


Reactor core diameter, По) ft . . ....:....... o 2.0 апа 2.5 
Volume concentration of fuel-element material, fg . . . . 0.06 and 0.12 
Fuel-element plate thickness, in. . . s ea a s e oea s. «© «© . 0.012 
NaOH coolant velocity, ?%%/вес................... 15 


Figure 15 is for МК ОҒ 100,000 pounds; figure 16 ів for Wr of 
150,000 pounds. Four separate plots are included in each figure, one for 
each of the four flight conditions. 


The NaOH coolant temperature leaving the reactor is included in the 
figures &s &n indication of the limit to which the maximum fuel-element 
temperature may be reduced. It is recalled that an assumed temperature 
difference of 100° F between the average NaOH temperature in the reactor 
and the effective wall temperature in the air heat exchanger has been 
incorporated in the data of figures 15 and 16. 


For a given reactor configuration and allowable maximum fuel-element 
temperature, figures 15 and 16 provide a value of reactor heat release 
required for any of the four flight conditions considered. Corresponding 
values of gross airplane weight may be obtained from ае 15 апа 14 
with Тұ, еғҒ as the linking variable. 


Figures 15 and 16 indicate the following general results for the 
two values of WK considered: 


(а) At subsonic flight (Mach number 0.9), maximum fuel-element 
temperature сап be maintained at 12009 Е and lower for reactor core 
diameters of 2 feet and volume concentrations of fuel-element material 
of 6 percent with reasonable values of reactor heat release and airplane 
gross weight. At 12000 Е, these values are tabulated_as follows: 
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Мк | Altitude Н We 
(15) (ft) (kw) (1b) 
100,000| 50,000 | 68,000 | 198,000 
100,000 | 30,000 | 59,000 | 160,000 
150,000 | 50,000 | 106,000 | 505,000 
150,000| 30,000 90,000 | 242,000 


Variation in reactor core diameter and volume concentration of fuel- 
element material has relatively little effect on maximum fuel-element 
temperature for a given reactor heat release. А maximum fuel-element 
temperature of about 1100° F can be attained by an increase in reactor 
core diameter to 2.5 feet or by an increase in volume of fuel-element 
material to 12 percent with operation at the following conditions: 


Ук Altitude H We 
(1b) (ft) (xw) (1b) 
100,000 50,000 72,000 205,000 
100,000 50,000 61,000 165,000 
150,000 50,000 110,000 510,000 
150,000 50,000 95,000 246,000 

For the subsonic flight conditions it is noted. that the increase in 
reactor heat release and airplane gross weight, resulting from а decrease 
in maximum fuel-element temperature from 12009 to 11009 F, is small. 


There appears to be а small advantage in reactor operation at uniform 
element temperature over operation at uniform power generation. 


(b) At supersonic flight (Mach number 1.5), maximum fuel-element 
temperatures must be of the order of 15009 to 17009 Е Рог а reactor core 
diameter of 2 feet and volume concentration of fuel-element material of 
6 percent. Increase in reactor core diameter to 2.5 feet and.volume of 
fuel-element material to 12 percent lowers fuel-element temperature to 
the 1200° to 1300° F level. For these conditions and a maximum fuel- 
element temperature of 1200° F, the following tabulation is made: 


Altitude H 
(ft) (kw) 
50,000 | 340,000 
30,000 | 260,000 


50,000 590,000 | 500,000 
50,000 410,000 | 275,000 
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Reference to figures 15 and 14 indicates that operation at the 
lower values of air heat-exchanger effective wall temperature seriously 
penalizes the airplane gross weight and required heat release for the 
supersonic 50,000-foot-altitude condition, but that flight at а Mach 
number of 1.5 at an altitude of 30,000 feet requires a larger heat 
release but a relatively smaller gross weight which is comparable with 
subsonic requirements. 


Summary of airplane conditions with low-temperature turbojet 


cycle. - For convenience in evaluation of the airplane engine and reactor 
requirements for the four flight conditions considered, a summary is pre- 
sented in table V of some representative results obtained in the airplane 
cycle performance, heat transfer, and reactor criticality study. 
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The. turbojet operating conditions and airplane component weights 
listed are based on engine operation at an air heat-exchanger effective 
wall temperature of 1000° F at all flight conditions. A difference of 
100° F between the average NaOH reactor coolant temperature and the air 
heat-exchanger effective wall temperature is maintained. 


The actual total engine air flow to the turbojet engines at each 
flight condition is listed; in addition, the total engine air flow 
corrected to static sea-level conditions is included. (The corrected - 
air flow is used to specify the &ir-handling capacity of & turbojet 
engine.) From designs of existing engines, it appears that a corrected 
air flow of about 300 pounds per second is obtainable with a turbojet 
engine approximately 4 feet in diameter. The corrected total engine air 
flow therefore serves as a measure of the number of turbojet engines 
required to maintain a given flight condition; these values are also 
listed in table V. 

Six and nine engines of the size specified are required to maintain Е 
flight at 50,000 feet altitude and Mach number 1.5 for respective values 
of Wr of 100,000 and 150,000 pounds. This compares reasonably well with 
the numbers of engines required for both subsonic conditions. Supersonic 
flight &t 50,000 feet altitude, however, requires on this basis 20 engines 
for a Wr of 100,000 pounds, апа 29 engines for а Мк of 150,000 pounds. 
The large number of engines required for this flight condition may impair 
aerodynamic characteristics and reduce obtainable lift-drag ratios below 
the assumed value of 6.5. 


For supersonic flight at 50,000 feet altitude, the air heat-exchanger 
effective wall temperature would have to be increased to about 1600° F, 
in order to reduce the required air flows to values comparable with the 
other flight conditions. 


Two reactor geometries of diameters 2.0 and 2.5 feet have been 
selected for operation at each flight condition. Both configurations 
provide for relatively low maximum fuel-element temperatures consistent 
with the attainment of reasonable airplane gross weights and reactor 
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heat releases. The uranium investments required Рог these reactors are 
below 50 pounds. | Пре AN | 


SUMMARY OF RESULTS 


The effects of reactor variables on criticality and maximum reactor 
fuel-element temperatures for sodium-hydroxide-cooled, moderated, and 
reflected reactors of length-diameter ratio 1 are presented for subsonic 
and supersonic propulsion of nuclear-powered aircraft. The following 
results, based on turbojet engine cycle operating conditions optimized to 
give minimum airplane gross weight, were obtained from the study for two 
values of weight of shield plus reactor plus pay-load plus auxiliaries 


(designated as Wk) representative of the divided-type shadow shield and 


the unit bulk shield, respectively: 


l. Flight at a Mach number of 0.9 at altitudes of 30,000 and 
90,000 feet for values of Wr of 100,000 and 150,000 pounds may be 
maintained with maximum reactor fuel-element temperatures of the order 
of 11009 to 12009 Е with а reactor core diameter of 2 feet. 


г. Ап airplane gross weight of about 200,000 pounds and reactor 
heat release of about 70,000 kilowatts are required for subsonic flight 
for Wr of 100,000 pounds. Gross weight and reactor heat release are 
proportionately increased for a value of Wr of 150,000 pounds. 


5. Flight at a Mach number of 1.5 at an altitude of 30,000 feet 
for both values of Wg may be maintained with maximum reactor fuel- 
element temperatures of the order of 12009 to 1300? Е with а reactor 
core diameter of 2 feet. At an altitude of 50,000 feet апа а Mach 
number of 1.5, similar maximum fuel-element temperatures may Бе maintained 
with а reactor core diameter of 2.5 feet. 


4. At 350,000 feet &ltitude, &irplane gross weight for supersonic 
flight is comparable with gross weight requirements for subsonic flight 
(&bout 180,000 pounds for Мк ОҒ 100,000 pounds); reactor heat releases 
required, however, are about four times greater than for the subsonic 
case (270,000 kilowatts for МК of 100,000 pounds). At 50,000 feet 
altitude, airplane gross weights and reactor heat releases required 
for supersonic flight are about 300,000 pounds and 350,000 kilowatts 
respectively, for Мк ОҒ 100,000 pounds. Gross weight and heat release 
are proportionately increased for Wr of 150,000 pounds. 


5. Enriched uranium investments for the hot unpoisoned reflected 
reactor, containing sufficient high-nickel-alloy fuel elements to provide 
appropriate heat-transfer surface for the aforementioned reactor and 
airplane flight conditions, are of the order of 35 and 50 pounds for 
core diameters of 2.0 and 2.5 feet, respectively. As fuel-element 
structural material in the reactor is reduced to zero concentration, 
the uranium investments approach 15 and 20 pounds for the 2.0- and 
2.5-foot core diameters, respectively. The excess uranium required to 
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counteract burnup апа the poisoning effects of equilibrium xenon, 
samarium, and other fission products, resulting from reactor operation 
for 24 hours at 300,000 kilowatts, is estimated to be less than 10 pounds. 
At thís poisoned condition, the temperature coefficient of reactivity for 
reactors with high-nickel-alloy fuel elements was calculated to be 
negative and of the order о? -0.00006 per ОЕ. 


lewis Flight Propulsion Laboratory 
National Advisory Committee for Aeronautics 
Cleveland, Ohio, July 28, 1952 
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APPENDIX А - SYMBOIS 


reactor coolant flow area 


-coolant specific heat 


neutron diffusion length 

reactor core diameter 

reduction in reactor diameter due to reflector 
equivalent diameter of coolant flow passage 
neutron energy 

airplane net thrust 


volume fraction of structural and fuel-element material in 
reactor core 


reactor heat release 

reactor heat-transfer coefficient 

neutron multiplication constant 

thermal conductivity of fuel-element material 
coolant thermal conductivity 

reactor core length 

mean square slowing-down distance for fast neutrons 
mean square diffusion distance for thermal neutrons 
airplane lift-drag ratio 

Maxwellian neutron density (per unit speed range) 
reactor atom density 

resonance escape probability 

neutron slowing-down density 

ratio of hydrogen to uranium 235 


reactor core radius 
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reactor radius 

reactor heat-transfer surface area 
fuel-element plate spacing 

temperature 

average reactor coolant temperature 

coolant temperature entering reactor 

coolant temperature leaving reactor. 

air heat-exchanger effective wall temperature 
maximum fuel-element temperature 

maximum fuel-element surface temperature 
plate-type fuel-element sandwich thickness 
mene Lor reflector pülecrness 
logarithmic neutron energy 
coolant velocity Е | 
neutron velocity 


airplane gross weight 


shield, reactor, pay-load, and auxiliary equipment weight 


airplane structural weight 

enriched uranium investment 

excess uranium requirement 
turbojet-engine air flow 
extrapolation distance 

macroscopic transport mean free path 


coolant viscosity 
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neutrons produced per fission 


average logarithmic energy loss per collision 


coefficient of reactivity 

absorption сговв section 

absorption cross section of structural material 
fission cross section 

Scattering cross section 

total сговв section 


Scattering cross section 


3 

coolant density 
яр temperature 
ZA macroscopic 
oath macroscopic 
Хт macroscopic 
Ls macroscopic 
Хт шасговсоріс 
ба microscopic 
Фф neutron flux 
pt adjoint function 
Subseripts: 
0 reactor core 
l reactor reflector 
m moderator 
th thermal neutron group 
f fast neutron group 


54 Sus МАСА RM Е52119 


APPENDIX B - TWO-GROUP DIFFUSION THEORY 


Two-group two-zone diffusion equations. - The neutrons are separated 
into two groups: а fast group consisting of all neutrons from fission 
energies to thermal energy, апа а thermal group consisting of all slowed- 
down neutrons in statistical equilibrium with the temperature of their 
surroundings. 


For the fast group, terms representing leakage, absorption, slowing 
out, and production of neutrons in the reactor core are given by 


^ 
Leakage: - 00 Аре 0 


Absorption: ЖА 7,0 $2,0 


Slowing out: 25 Р 9 Pe,o 
Production: Ўр, 69 с + Хр, tho УФ+ћ,0 


For the thermal group, terms representing leakage, absorption, and 
production of neutrons in the reactor core are given by 


е TR th,.O 
• --- ‚О 


Absorption: sas Pih,o 
Production: %5,%,0 Фе, о 


For the reflector, terms representing leakage, absorption, and 
slowing out of the fast group are given by 


[3l 
Leakage: --3- A Pe, 


Absorption: Zar, Фет 
Slowing out: Ба, Р,1Фе,1 


Terms representing leakage, absorption, and production of neutrons 
for the thermal group of the reflector are given by 


ATR, th,1 

Leakage: - a zc m AQ th,1 

b : 2, 

Absorption ША, th,19tn,1 

Production: às 0,19 7,1 

The macroscopic cross sections and mean free paths appearing in the 

foregoing terms must be appropriately averaged over the energy range 
constituting the pertinent neutron group so that the various terms 


represent the correct distribution of neutron processes occurring within 
the reactor. 
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For each group and reactor zone, expressions for the conservation 
of neutrons may be written 


- Mm £20 ^ Фес + (Xa, e,o + 25,2,0) 9,0 - “в,ғ,0 99-07 
Хғ,ыҺ,09%7ьҺ,0 = © (21) 
и Amy th,0 A^ 94,0 + Кац Ф+ъ,о - 28,2,0 ?ғ,0 = 0 (B2) 
CM АФ», у + (Ёд, Р, + 25,2,1) 95,17 0 (83) 
= mama AGin,i + Za, th, 19 tn,1 7 25,ғ,1%9ғ,1 © (ва) 


These equations тау be rearranged as follows: 


5(% де Ж, ) У D 
А,Ғ,0 * “s,£,0 5 F,£,0 
App oer c SS We Pe ot 
í Атв,?,0 £,0 Xm,r,o 2А,2,0 01:9 1 $59 
5  ÀTR,th,O 2F,th,o” 
Хтв,ъһ,о Әтв,Ғ,0 ~A,th,o Аз 20 148,0 у 
^ _ SA th,0 ф " 5%5,Ғ,0 NTR,f,0 Peo = 0 (B6) 
th,0 7 km, tn,0 00 7 Am,r,o ATR, th,o +? 
5(2д.р,1 + 28,2,1) 
А Фе, - SMS CORR MG E 2-2 Op 1 = 0 (B7) 
3+ 3 
5% 52 X 
A, th, S,f,1 ÀTR,f,1 
^ p - 2 2 Ф + 2 „эел. ф = (B8) 
бб, ,th,l св, 1 ATR 2,1 TR th,l = 


The mean-square в оф пр Жаы» Length Le 0,1 and mean-square 
thermal diffusion length L°tnh,o0,1 are defined as 


|2 _ _ Атв,ғ,0,1 
1793 5 ҰАР,0,1725,6,0,1) 
^ 
L = LL PL 
өшін ЗА th,0,1 
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Inasmuch as the macroscopic cross sections are the effective averages 
over the pertinent energy group, the resonance escape probability may be 
suitably approximated in terms of these effective values. 


: В 55,ғ,0,1 
5Һ,0,1 A,f,0,1 + 25,2,0,1 


15 follows then that 
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Мтв,Ғ,0,1 ‚ (® £051: Ус £,0,1 1e 0,1 


; ТР 
> ni tos губу кре == 
S; T;0,1 S f,0.1 di Pth,0,1 


2 
Nm,f,0,.1 (Фл,ғ,о,1 * 25,2,0,1\ _2 U£,0,1 
= L"£,0,1 


? 2 = 
ЛР. Та род (1-Ptn,0,1) 


The multiplication constants Кр and Kt, representing, 
respectively, the number of neutrons born per и absorbed in each 
energy group are defined as 


к,г о, 
Kp ХА, г, x 
ка = ФЕ th Os 

A,th,O 


With these о. the зрее equations may be rewritten. 
For the core, 


1-Ke(1-Pq& 0) А tho Еһ 
А Фе с 7 ха > Фғ,0 EE s Тат Ф+һ,о = О (B9) 
Т, f,O Тн,І,0 L $6,0 
АФ — qun,o + уина. 20 0 0 (210) 
th,O0 - 72——— 9,0 -5 — Peo = 
, Leino | TR,th,O L^. 0 4 


For the reflector С 
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à PII 2 Фе,1 О (B11) 
2 
А р 
1 TR,f,1 Pth,1 
Аф - —5——— Ф 2-2 —$ = 0 B12 
th,1 ы. thl + Kop thi тё tt (B12) 


Boundary conditions and solutions. - The preceding simultaneous 
linear differential equations are solved subject to the following 
boundary conditions for the case of spherical geometry. 


GeSe 


The continuity of neutron current at the core-reflector interface 
requires that 


(B13) 


AR th,1 
BUE acia и = =. eee 
2 юра, 5 УТ, ве 


The continuity of neutron flux at the interface also requires that 


Maud 
ВО ева te тан, 


(814) 
Pth, 1 ма. Фън, || 


At the outer face of the reflector, the neutron fluxes must go to 
zero at the extrapolated boundary; it is assumed that the extrapolation 
distance is the same for both fast and thermal neutrons and small in 
comparison with the reflector thickness so that 


= О 
j f | ге Набу 


(В15) 
T E 


r=Rettr 


The neutron fluxes must be finite everywhere in the reactor, 
that is 
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Фе Ей у < œ 


(816) 
Penlan г < 9 


The general solutions of the pairs of simultaneous equations (9) to 
(12) are obtained by letting Аф= айф where ой ів the same constant 
for the fast and thermal fluxes in each zone but differs for the core and 
reflector zones; a? is designated as ag for the core апа ав а Рог 
the reflector. Introduction of АФо,1 = ao, 190,1 in equations (9) to 
(12) yields quadratic expressions for ао and a4? in terms of the core 
and reflector compositions, respectiveiy. The general solution for the 
fast or thermal flux in each zone is then a linear combination of the 
general solutions of АФ = ap corresponding to the two roots of the 
quadratic appropriate for each zone. 


Application of the foregoing boundary conditions to the general 
solutions for the fast and thermal flux in each zone leads to four homo- 
geneous linear equations involving four arbitrary constants and the 
unknown critical radius of the reactor соке апа thickness of the reflector. 
The vanishing of the determinant of coefficients of these four equations 
for an assigned thickness of reflector is the condition for the evaluation 
of the critical radius of the core. 


With the critical dimensions of the reactor known, three of the 
arbitrary constants are solved for in terms of the fourth. The flux 
distributions in both the core and reflector are therefore determined 
except Рог ап arbitrary multiplicative constant. Detailed mechanics of 
the two-group method аге presented in references 17 апа 18. 


The. reflected reactor sizes calculated for spherical geometry аге 
translated to cylindrical geometry in the following manner: 


(a) For the unreflected reactor, the critical sizes are determined 
by the relations 


2 
«с^ = ла for spherical geometry 
(В! +6) 
2 nê 2.4052 
+ LL for cylindrical geometry 


ТО = (наб) (Ras)? 
where the buckling constant о із а function only of the composition 
of the unreflected reactor and hence is independent of the geometry of. 
the system. (Но is here the cylinder radius, R'e the sphere radius, 
and 5 the extrapolation distance at which the flux becomes zero.) 
Hence, for the same reactor composition, 


2020 
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ле » де " 2.4052 
(В! +6)2 (Н+25)2 (Вс+5)2 


or, if the cylindrical reactor is of length-diemeter ratio equal to 
unity He = 2Re, 


Re + 5 = 0.9144 (Н! +5) 


B ~ 0.9144 


(b) For the reflected reactor, the assumption is made that the 
reflector savings (difference in core radii for unreflected апа reflected 
reactors) is the same for both the spherical and cylindrical geometries. 
Hence the reflector savings calculated for any given core and reflector 
thickness for the spherical geometry can be directly applied to the 
unreflected cylindrical critical dimensions to obtain the reflected 
cylindrical critical dimensions. 


Effect of voids on criticality of bare reactor. - The introduction 
of uniformly distributed void space in a bare reactor reduces the densities 
of all constituents in the reactor by the factor (1-f.) where fy, is 
here the fraction of void volume. Hence, all macroscopic cross sections 
&re reduced by the factor (1- fy). The relation between agé and the 
reactor nuclear constants (see references 17 and 18) shows that о is 
then reduced by the factor (1-Р,) for the same relative core composition. 
Hence, the critical dimensions of the bare reactor are increased by the 


factor l that ів. 


CN)! 
D [wath “ода 1 


D [no void] = (1-#.) 


For the same relative core composition, the uranium investment 15 
directly proportional to its reactor concentration and to the reactor 
volume; that is, 


p? [with void| 1 
| x 2 
р? [no void] (1-f.) 


WU [with void] 
W no тоза 


Effect of voids on criticality of reflected reactor. - The assumption 
is made that the reflector savings are unaffected by introduction of small 
amounts of void volume in the core ОҒ а reflected reactor. 
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APPENDIX C - EVALUATION OF TWO-GROUP THEORY CONSTANTS 


The constants required for solution of the two-group equations are 
as follows: 


ne mean-square slowing-down length 
i^ mean-square thermal diffusion length 
Ата, fast-transport mean free path 


hor, th thermal-transport mean free path 


Pth resonance езсаре probability 
Ke fast multiplication constant- 
Ken thermal multiplication constant 


Ihese constants &re separately evaluated for each core and reflector 
composition. Multiplication constants are, of course, zero for the 
reflector. | | 


Ihe fast parameters were averaged over the energy spectrum of 
neutrons born during fission. The thermal fission neutron spectrum 
was taken from reference 5 and normalized to unity over an energy range 
extending from 25,000 ev to lO Mev. 


Total eross sections for Na, O, В, Ni, апа Fe ав а function of 
neutron energy were taken from reference 8. Resonances for these nuclei 
occur only at 5000 ev and above; for these energies the neutron scattering 
width is very much larger than the neutron absorption width so that these 
resonances have been taken as substantially scattering. The values used 
for total cross section for H were those obtained from measurements in 


H50. 


Thermal absorption cross sections were taken from reference 19; 
the pile oscillator values were used wherever available. The variation 
Y absorption cross section with neutron energy was assumed to follow the 
l/v law. пещ | | | 


Fission, absorption, and scattering cross sections for K-25 end 
product uranium mixture (91.5 percent 0255, 1.5 percent U254, 7 percent 
0258) normalized per atom of 1239 were taken from references 20 and 21. 
Metal densities were taken from reference 22. A summary of cross sections 
and densities used in these calculations has been presented in table I. 


2525 


себе 


КАСА RM #52119 HERS 41 


Fast Parameters 


Mean-square slowing-down length, ІЗ». - Ав previously mentioned, 


mean-square slowing-down lengths were calculated for hydrogenous mixtures 
by the method of reference 6 іп which & formula, equation (108), for 


та, is derived taking into account the scattering but not the slowing- 
down properties of elements other than hydrogen. This introduces small 
error provided the macroscopic scattering cross section of hydrogen is 
greater than the macroscopic scattering cross section of the heavy 
elements in the medium. A limiting form of the rigorous formula, 
equation (A3) of reference 6, was used for the present calculations 
Since it provides results in good agreement with the more rigorous 
equation (108) and involves considerably less computational labor. 


The values of 12. for water of unit density, by equations (108) 


and (А5) of reference 6 averaged over the fission spectrum, were calcu- 
lated ав 26.9 and 25.4 square centimeters, с е; These values 
compare with the experimentally determined value of р for water of 

55 square centimeters of reference 25. It is customary in reactor 
calculations to correct the caiculated value of тёр to the experimental 


value. This multiplicative correction for water from the I^. ав caleu- 
l&ted by equation (А5) of reference 6 is 55/25.4 = 1.30. 


Inasmuch as no experimental data are available for NaOH, the 
correction factor to the calculated value of 2, for Маон was determined 


by ascertaining the mixture of H20. and №20 corresponding to NaOH. 
It was found that а mixture by volume of 40 percent H20 of specific 
gravity 1.00 and 60 percent На?0 о? specific gravity 2.27 corresponded 
to NaOH of specific gravity 1.77. Assuming the correction factor to be 
proportional to the percentage of Н20 in the mixture, a correction 
factor of 1.12 was established to apply to all values of Ip as 
computed from equation (А5) of reference 6. These values of Тер were 
rther corrected to the operating temperature of interest by assuming 
L°e to be inversely proportional to the square of the density of NaOH. 


Macroscopic cross sections. - The remaining fast parameters were 
calculated from the energy distribution of neutron flux, given by age 
theory for an infinite medium of the same composition, as the weighting 
factors on the energy-dependent cross sections (see references 18 and 24). 
For example, the effective fast macroscopic absorption cross section 
2a,f is given by 


> (а) Ф(0) du 
> Uth 


А,? 
| p (а) du 
"th 


(C1) 
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where the flux $'(u) is obtained as 


q(u) 
ф'(ц) а EZ.) | (c2) 


and where q(u) is the age-theory neutron slowing-down density through 
any energy interval in an infinite medium. 


Similarly, the effective fast macroscopic fission cross section 
2, and transport mean free path ATR, ғ may be evaluated as 


0 
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2 'p(u) ф' (u) du 
Чън 


Lor - б (C3) 
| ф'(0) du 
"th | | : 


О 
| Уан" (a) Ф' (а) аш 
wth И 


TR,f " О 
Ф'(м) du 
Uth 


A (C4) 


The fast multiplication constant is therefore given by 
Ке = vip, г/ де: 


The resonance escape probability py, ів the value of q at 
thermal energy мр; inasmuch as а fission spectrum normalized to unity 
has been employed. | - 


The integrals required аге evaluated numerically over the entire 
fast energy region and in suitably smail energy intervals. 


Thermal Parameters 


Macroscopic cross sections. - The thermal parameters of both core 
and reflector involve macroscopic cross sections which must be suitably 
averaged over the Maxwellian distribution of neutrons in the thermal 
region to represent their effective values. For example, in order that 
the term БА, th?th represent the true rate of neutron absorptions 


occurring per unit volume of reactor, BA, th must be some effective 
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average value if Фи is given ав the thermal flux of neutrons corres- 
ponding to the most probable speed of the Maxwellian distribution 
corresponding to energy Eth. 


Ihe effective value of 2a,th averaged over the Maxwellian distri- 
bution of neutron flux  vM(v) Тв given by 


Фо 
2, vM(v) dv 


А | 
ccn don uteri (c5) 


vM(v) àv 
О 


For absorption cross sections which follow the 1/у law, ав is 
usually the case for thermal neutrons, the value of È A,th evaluated 
from equation (C5) turns out to be exactly the local value of Од 


corresponding to the average speed of the Maxwellian distribution of 


neutron density M(v). Inasmuch as the average speed is 2 = 1.128 
x 

times greater than the most probable speed, the effective value of 

La tn is 0.886 times the value at the most probable energy to which 

it 1s customary to refer measured cross sections. (The pile oscillator 

values of thermal-absorption cross section used herein have all been 

referred to the accepted value for gold at the most probable neutron 

thermal energy). 


The same results apply to the macroscopic fission cross section of 
uranium Хр, which very closely follows а l/v variation for the 
thermal region. 


Similar considerations are made in obtaining an effective value 
of thermal-transport mean free path ATR, th: The effective value of 


АТВ th is given by т 


Mop vM(v) ах 
О 


ATR, th = де 
vM(v) àv 


(c6) 


O 
where Лов! is the local value for neutrons of а particular velocity. 
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For any particular neutron velocity, Z', and  A'pg are evaluated 
ав the sum of the products of the reactor atom density of each constituent 
эпа the pertinent cross section: 


T 
LA = Мі 6д,1 + No бд,2+...+ Na бА, а (ст) 
1 
NOE = Ny Ств,1 + No бтр,2 t+... + № OTR,n (са) 


Effect of chemical binding. - The effect of chemical binding of 
the hydrogen in the molecule for both water and NaOH alters the angular 
distribution of neutrons scattered by hydrogen. Since this effect - 
occurs almost exclusively in the thermal region and varies rapidly with 
neutron energy, the average value of the cosine of the scattering angle 


cos 6 ав a function of local neutron energy must be known. 


In the absence of these specific data, Radkowsky in reference 7 has 
made use of the "Born approximation" (reference 25) in evaluating an 
effective value of. А, the atomic mass of the scattering atom, for use 
in the formula for сов 6 x 2/34 applicable for isotropic elastic 
scattering in the center of mass system. The "Born approximation” 
indicates the scattering cross section to be proportional to the square 
of the reduced mass of the neutron and scattering atom. For neutrons of 
energies well above the region of chemical binding, the hydrogen atom 
is effectively free and presents a mass of unity and os of 20 barns. 
Hence, the following proportionality may be formed from which an 
effective value of A 148 determined: 


A y 
С 1 — 
5 Ari 
20- = "i74 (c9) 


where ба" is the local value of hydrogen scattering cross section 
Observed in measurements on water. 


_____Ву means of equations (C6), (C8), (C9), and the formula 
сов Ө = 2/5А, Radkowsky has checked the experimental values of- thermal 
diffusion length for water for a range of temperature. Effective values 
of Аъ + for the hydroxide reactors herein considered have been 
numerically evaluated by the same procedure. 


The remaining thermal parameters are evaluated as follows: 


The thermal multiplication constant Ку ів given by 


2525 
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The mean-square thermal diffusion length is given by 


га, = 1 NTR, th 
В 3 Esth 
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APPENDIX D - TEMPERATURE COEFFICIENT ОҒ REACTIVITY AND EXCESS 
URANIUM REQUIREMENTS 


The derivation of the first-order perturbation formula presented in 
appendix I of reference 12 for the two-group two-zone problem has been 
modified to include the fast fission effect. This modification involves 
no changes in the derivation procedure of reference 12 and, for this 
reason, the details of the modified derivation sre omitted herein. 


The perturbation formula gives the effect of small changes in the 
nuclear properties of any portion of the reactor volume on the reactivity 
of the entire reactor. The perturbation formula is herein applied to the 
problems of determining reactor temperature coefficient of reactivity and 


excess uranium requirement to overcome burnup and fission-product 
poisoning. 


Perturbation formula. - The two-group neutron diffusion equations 
for the reactore core are 


редфе - (D, & + 2g ¢) Фр + Ки Баљу Ptn + Ke Za, e Pp = 0 (D1) 


ПаһАФьһ ~ ZA,tn Pth + Pg,e Фр = 0 (D2) 


where the subscript 0, used to indicate reactor core, has been dropped 
in this section. 


In matrix form, these equations are represented by Hp = 0, where 
the matrix operator Н is defined as 


БА - 24 tn Жа,ғ 


Я = 
Кер ФА, th рел - (Ед,р + 25,2) + Ke дог 


and Ф is the column matrix 


Фън 


Фе 


The operator Н defined herein includes the fast fission term 
Ke 5д,ғ which was taken as zero in reference 12. 


2525 
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The adjoint matrix H*, obtained by reflection of &11 elements of 
H about its principal diagonal, is given by 


БА - 2 th Кен ZA th 
H* = 


У, ГА - (X, p + Бар) + Kp Бр, 


The adjoint matrix equation defining the adjoint functions is then 
H*p* = О and the adjoint equations are given by 


АФ» - (Z, + Bsr) ф + Ke ba gU + е? = 0 (D3) 


ТА th - Za tn P'en + Kin Ben = 0 (D4) 


Ihe addition of the term Kp ФА, е in the second-row second-column 
element of the operator matrices Н апа H* introduces no change in the 
detailed derivation of the perturbation formula given in appendix I of 
reference 12. Hence, only the final formulas &re presented herein. 


The perturbation formula, written in terms of the neutron fluxes 
Фън and Фр and their adjoint functions Pih and tp, is: 


5, Spl | € tn 
(9 tn 9 г) ат 
reactor бе бај | e 
P шы ышы M — — = 7. 


(Kin ZA tn Фън Pop + Ке дог Фе 9^ ejat 
core 


where 


С? 
| 


а = 5 Е -2а, tnl 
ъ= 5 LO 
бе = 5 [Ken “А, en] 


ба ж б [De - (Za г T 25,7) + Кр SI 


SEND 
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If the indicated matrix operations are performed and P expressed 
in terms of the various contributions to the reactivity due to the 
individual changes in the pertinent nuclear constants, the following 


results аге obtained: 


Р = P(Dtn) + Р(ре) + P(A, к) + P(Zs р) + P(Z, tn) + РР e) + Р(РЕр th) 


where 


Р(Бан) = 


Р(р2) = 


(2А ғ) = 


P(g ғ) 


P(ZA, th) 


P(VZp, th) 


P(v2y g) 


(D6) 
- | VO" ыы "УФ 5(Dgs)àt 
vol (D7) 
Y 
+ 
- | УФ p 7 У%;, 5(0.)ат 
vol y (D8) 
= P tp ФР 5(2 pjat 
vol 
D9 
= (29) 
| $ Фе - Ф”ьФ,)5(2а јат 
УОТ. Я | (510) 
Б Pth 5, )ат 
vol 
= (011) 
Фе Pin Әр, въ)л 
vol REN (р12) 
Y 
Ф.Ф, S(v Ep e)àT 
и (D13) 
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+ 
Ү = (ев Фъ Фе + 92р, 7:9 „јат 
соге 


Adjoint equations. - For use in the perturbation formula, the 
adjoint functions ФТр and Prep must be known as a function of 
position (radius r, in this case). 


Equations (3) and (4) are the core adjoint equations. 
= + + 
De oA 2,0 Č л,г,о + 28,2,0) Ф £50 + Кр 5ААР,О $ + О 2а,р,о Pth, = 0 
(D3) 


+ + + 
Dtn,o^Q tho - Pa,tn,0 9 tho + Étn д, tho Ф ғ,о = © (D4) 


The reflector adjoint equations, obtained in the same manner as 
for the core, are 


+ + + 
De, 49" е, - (E e, + 25,2,1) V 2,1 + 23,217 th = 0 (DIA) 


+ + 
Den, AP th,1 - 4A, tn Ф th, = O (015) 


In the derivation of the perturbation formula, the same boundary 


conditions are applied to the adjoint functions as are applied to the 
fluxes, namely: 


+ + 
Ф ғ,0= Ф £,1 
+ + 
Ф tho =? thi 
at core-reflector 
De 0792,0 = De, N? 41 interface 


+ ua 
Din, 079 th,o = Ptn, V? th,1 


Фет = Q'uni- О а reflector outer surface 
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t + 
М 3 ipu ци r 


The solution of the adjoint equations subject to the foregoing 
boundary conditions is identical, in procedure, to the solution of the 
neutron flux equations subject to the same boundary conditions. 


Temperature coefficlent of reactivity. - The procedure used for 
determining the differential changes in the pertinent nuclear parameters 
with temperature, required for use in the perturbation formula, is pre- 
sented herein. 


2525 


For the fast parameters, two principal effects of change іп temper- 
ature аге accounted for: 


Effect A - accounts for change in neutron energy range constituting 
the fast group. | | | 


Effect В - accounts for change in density and hence in density of 
nuclei of NaOH in reactor core only inasmuch ав uranium and structural 
material remain fixed in the reactor. 


The relations for the fast parameters are 


d£, p dda Win 


ат dut, ат (D16) 

dp, Ёр, йш (ті) 
ат T duth ат 

dàm, f Мн? Mth А dom (D18) 

ат а ат 8,1 сат 

ады р don 
2 [o 

2—6 Уве р-ат (D19) 


In relations (016) to (D18), the respective values of à/àucg are 
obtained from equations (Ci), (C3), and (C4); the subscript m refers 
to the moderator. 


Effect B is neglected in equations (D16) and (D17) because the 
fast absorption by NaOH is small. Effect. A is negligible in comparison 
to effect B in the neutron slowing-down process and hence is neglected 


себе 
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in equation (019). Effect B is included in equations (D18) апа (D19) 
in а manner implying that the diffusion апа glowing-down processes are 
due primarily to the Маон. 


For the thermal parameters, two principal effects of change in 
temperature are accounted for: 


Effect B - same as for fast parameters. 


Effect C - accounts for change in microscopic cross sections with 
change in thermal temperature. 


For all nuclei within the reactor excepting Хе155, the microscopic 
absorption cross section сд is assumed to follow the l/v Лам in the 
vicinity of thermal energy. Hence, for these nuclei бд is inversely 
proportional to AT во that 


dcA th ат 


TA, th er 


For Хе155, в quantity сХе which is а function of temperature is 


defined аз 


„Хе а сдХе 


The relations for the thermal parameters аге 


On th арт асты th 


dX. th "m 1. Хе „Хе 1 
— Абы e 
= Za th,m pat + Za, th (= + № (021) 
24654 ХА th) l 
2 = — 
ат = Ken Фа ъп =) (р22) 


Тһе foregoing changes in nuclear parameters аге assumed to be 
homogeneously distributed over the entire reactor core. Мо changes 
in the reflector parameters are taken as contributing to change in 
reactivity. | 
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Estimate of excess uranium requirements. - Use of the perturbation 
formula to calculate the additional uranium required to counteract the 
effects of fission-product poisoning апа uranium burnup is described 
herein. The principal poisons are Хе155 and 8ш149; a host of other 
fission-fragment poisons appear &nd &re accounted for by the &ssumption 
that for every uranium atom fissioned the equivalent of а single poison 
atom appears with а microscopic absorption cross section ОРА th equal 
to 100 barns. у 


2525 


The excess uranium requirement is given by the relation stating 
that the reactivity decrease due to accumulation of fission-product poison 
and due to fuel burnup must equal the reactivity increase brought about by 
addition of extra fuel. 


The assumptions used in the calculation of these reactivity changes 
are as follows: 


(1) Introduction of the fission poisons affects only БА th: 


(2) Only the uranium thermal absorption and production parameters 
are affected by subtraction or introduction of uranium into the reactor 
core. This assumption implies that the reactor is principally thermal. 


(3) The fission poisons and added uranium are assumed to be distri- 
buted uniformly over the volume of the reactor core. 


Let: 


anv number of uranium atoms per cubic centimeter of reactor core 
burned up during the required reactor operation. 


WW uranium concentration required for the hot unpoisoned reactor 
(Investment corresponding to this concentration is given by 
the criticality calculations.) 


B fractional increase in uranium concentration, over that calculated 
for the hot unpoisoned reactor, to counteract reactivity effects 
associated with fission product poisoning and fue] burnup. 


The weighting factors in the expression for the change in 
reactivity P, are evaluated as follows: 


(а) To calculate the reactivity decrease due to poison, P@ л th) 
(equation (D11)) is evaluated by taking : 


е „Хе т 49m 
Sn tn) = Бес A,th * yan а A,th * (59) о? Ath 
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(b) To calculate the reactivity decrease due to fuel burnup, 
Р(5д,ыһ) (equation (D11)) is evaluated by taking 


B(Za th) = - (ӘМ) o7, th 
and РоЖ,вн (equation (D12)) is evaluated by taking 


(vie th) = - (ЁН) cp th 


(<) To calculate the reactivity increase due to the addition of 
uranium in excess of that required for criticality of the hot unpoisoned 
reactor, P(Z A,th) (equation (D11)) is evaluated by taking 


5(2л, th) = BNU oU, th 


and Ро р, th) (equation (D12)) is evaluated by taking 


S Ep tn) = vpNU OUP th 
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TABLE I - NUCLEAR PROPERTIES OF REACTOR MATERTALS 


НАСА 7" 


0.0258 X 1024 


Atom Thermal — M pem cM 
pem cM 
жеш = 


“Thermal energy Eth = 0.092 еу cor- 


responding to an average temperature 
of 14509 F. 


3696 


t 
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TABLE II - TWO-GROUP TREORY REACTOR CONSTANTS 


(а) Core. 


Reactor|R - Em, f,0 мањој 20 Е ав, Кен | Pth,o 
( ош) (em) |(сш@)| (em?) 
im 
200 
50 | 5.858 | 1.784 |105.6| 5.066|1.658|1.850| 0.555 
150 | 5.884 | 1.795 105.6|12.204|1.254|1.461| .776 
5.895 | 1.798 [105.6|18.769|1.104|1.121| .854 
га 4.087 | 1.902 |104.9| 9.54 11.665|1.852| 0.729 
4.125 | 1.910 1|104.9|28.765|1.285|1.552| .902 
5.766 | 1.723 | 89.2] 4.299|1.75011.961 0.549 
5.859 | 1.740 | e&9.2|27.611]i.215|1.247| .917 


(b) Reflector. 


Ма, TR, th,1 L^ и г th,1| Pth,1 
EU (ст) (em?) 


ГТ 
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TABLE III - NUCLEAR PROPERTIES OF FUEL-ELEMENT ALLOYING MATERIALS 


Effective? 
cross Sections cross sections thermal 


(0.025 ev) Y = Nc macroscopic 
0. cross sections 
(0-05 ег) (0.092 ev) 


LL 
6 
1 
4 
7 
7 
5 
5 


Әл тегадей over Maxwellian distribution of neutron flux. 
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TABLE ТУ ~ EXCESS URANIUM REQUIREMENTS AND TEMPERATURE COEFFICIENT OF REACTIVITY FOR  NaOH-COOLED AND 


AAT 


MODERATED REACTORS. OPERATING POWER, 500,000 KILOWATTS; OPERATING TIME 24 HOURS 


Core composition: NaOH, 0.82; 
Ni, 0.08; На, 0.10. 


| Reactor 111. Core composition: NaOH, 0.90; 
: Ма, 0.10. 


R - 400 
D, = 5.10 ft 


R = 100 
De = 1.52 ft 


tp = 6 in. 


W « 14 Ib 
AU ы 4.5 1b 


+.0ҳ2737 


~ .Ox2461 
= 04207 


+.04265 +.0 44 


- 0,191 


х ПИВО КСО А НИ 22 ИЕСИ 


Lath %.051691 мі іші ~ 04154 -.04194 | +. um -.05728 | 4.0.2591 
without 
сілеп | 


-.05116 -.0,526 
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TABLE Y - SUMMARY OF REPRESENTATIVE ROCLEAR-POWERED-AIHELAME DATA FOR VARIOUS FLIGHT CONDITIONS 


Turbojet operating conditions 
Air heat-exchanger wall temperature, Op 
Optimum turbine-inlet temperature, ОР 
Optimum compressor pressure ratio 


Reactor, shield, payload, auxiliaries, lb 
Engines ani ait heat exchangers, lb 
Airplane structure, 1b 
Airplane gross weight, lb 

Total engine eir flow, lb/sec 

Total engine air flow corrected to static 


sea-level conditions, lb/sec 
Number engines of 300 lb/sec corrected 


Fuel element sandwich thickness, in. 

NaOH coolant velocity, ft/sec 

Reactor length and diameter, ft 

Fuel-element volume fraction 

Equivalent diameter of flow ревазде, in. 
Fuel-element maximm temperature, °F 

NaOH temperature into reactor, Uy 

Nai temperature out of reactor, OF 

Uranium investment ко fuel elements), 1b 
Uranium investment (stainless steel elements), 15 


TEETER 
ac B EER 


аНевасіста are reflected by 6 inches of МАСИ and operate at uniform wall temperature. Difference between average 
heat-exchanger effective wall temperature taken as 100? F. 


:88 


г. 
e 
ез 
M 

сл 


МАНЕ 
ТТ 


5256 


99 


6т125 ИЯ YOVN 
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NaOH reflector 


М. 


NaOH Plate-type sandwich 
fuel elementis 


“УАУ А РУ 
| ARAN г: 
i 


ИИЯИЯ 
МАКА“ 


Figure 1. ~ Schematic diagram of NaO0H-cooled, moderated, and reflected reactor and typical 
heat-transfer configuration. 
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Shield 


NaOH reflector 
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NaOH pump 


Primary 
NaOH -to-1iquid- 
metal heat exchanger 


Shield 


Liquid-metal pump 


Exhaust 
diffuser compressor Liquid-metal- turbine nozzle 


Figure 2. - Schematic diagram of nuclear-powered turbojet engine. Primary coolant, 
NaOH; secondary coolant, liquid metal. 
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Ty ex (uniform heat generation) 
Ta ex (uniform wall temperature) 


с,ех 


erroow Wu VoVN 


Average NaOH temperature Те ay 
7 


Е = 
| Assumption: 


о 
Tear М Ту ег? SOME 


REN] 
P нш well temperature Ty err 
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Figure 3. - Schematic diagram of heat-exchanger system and relative temperature distributions. e 
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Enriched uranium investment, WU, 1b 


180 


 [LLLLLLLLLLLLLLLLLL. 
ШЫ 


Reflector 
eee | 


= of |X со LT Wy ||| 
ТОЛИ РАЦ 


Ratio of hydrogen 
to READER 


100 


LLL LL МО 
к yt (ДА ЖЖ дааа 
ТУК УДА || | ||| | | 
А77 11 бан Шш 
ш ЕЕ ERN 


Reactor core diameter, D 


ft 
с? 
(а) Reactor I. Core composition: NaOH, 0.82; Ni, 0.08; Na, 0.10. 


Figure 4. - Enriched uranium investment for Na&OH-coaled, moderated, and reflected Sri oe Length- 
diameter ratio of cylindrical core, 1.0; average moderator temperature, 14500 
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Епгісһей uranium investment, мо, lb 


‚8 1.2 1.6 2.0 2.4 2.8 5.2 5.6 4.0 4.4 
Reactor core diameter, D,, ft 


(b) Reactor II. Core composition: Маон, 0.82; Fe, 0.08; Na, 0.10. 


Figure 4. - Continued. Enriched uranium investment for NaOH-cooied, moderated, and reflected 
reactors, Length-diameter ratio of cylindrical core, 1.0; average moderator temperature, 
14509 F. У i 
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thickness 


= 


Enriched uranium investment, wi, lb 


20 


КИ Satta te кЕЕШЕ 
|І Не | | | | | || | | | || 
| || || А 
Reactor core dise По, ft 


(c) Reactor III. Core composition: NaOH, 0.90; Na, 0.10. 


Figure 4. - Continued. Enriched uranium investment for NaOH-cooled, moderated, and reflec- 
ted reactors. Length-diameter ratio of cylindrical core, 1.0; average moderator temper- 
ature, 1450° F. 
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Enriched uranium investment, WU, 1b 


ft 
(d) Reactor IV. Core composition: NaOH, 1.00. 


Reactor core diameter, Do» 


Figure 4. - Concluded. Enriched uranium investment for NaOH-cooled, moderated, and reflec- 
ted reactors. Length-diameter ratio of cylindrical core, 1.0; averege moderator temper- 
ature, 14509 F. 
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Наш 
Line of reflector savings 
equal to reflector thickness 
Ratio o drogen 
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Figure 5. -.Representative reflector savings for NaOH-cooled, moderated, and 
reflected reactors. Average reflector temperature, 1400? F. 
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Reactor 
composition 
by volume 


Маон 0.82 


Ni 
Na 


| TAT 
AVA е = ГГ 
“ТТТ LII 
И к=з 


Beal. Ieee 


Enriched uranium investment, м7, lb 


Reactor. core diameter, Па, ft 


Figure 8. - Comparison of critical reactors I to IV. Reflector thickness, 6 inches of NaOH; 
length-diameter ratio of cylindrical core, 1.0; average moderator temperature, 14509 PF, 
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Figure 7. - Generalized critical size and enriched uranium investment as a function of thermal neutron absorption parameter of reactor structure 
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Figure 10. - Radial distribution of fraction of fissions produced by thermal neutrons. 
14509 F; ratio of hydrogen to uranium atoms В, 100. 
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Figure 15. - Variation of required reactor heat release with air heat-exchanger 
effective wall temperature for optimum turbojet cycle. Data based on refer- 
ence 14. Ratio of structural to gross weight W,/W,, 0.50. 
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